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Abstract 

Charged current scattering of i/„ on 12 C has been studied using a ir + decay- 
in-flight Vn beam at the Los Alamos Neutron Science Center. A sample of 
66.9=L9.1 events satisfying criteria for the exclusive reaction 12 C(f M , N g , s , 
was obtained using a large liquid scintillator neutrino detector. The observed 
flux-averaged cross section (5.6 ± 0.8 ± 1.0) x 10 1 cm 2 agrees well with reliable 
theoretical expectations. A measurement was also obtained for the inclusive 
cross section to all accessible 12 N states C(y^ [i~)X. This flux-averaged 
cross section is (10.6 ± 0.3 ± 1.8) x 10~ 40 cm 2 which is lower than present 
theoretical calculations. 
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1 Introduction 



Neutrino-nucleus cross sections are needed for the interpretation of measurements by 
many low energy neutrino experiments as well as for modeling various astrophysical 
processes such as supernova explosions. Low-energy neutrino-nucleus cross sections 
are also of interest because of their role in nuclear structure studies. The cross sections 
contain contributions from both axial vector and polar vector nuclear currents and 
thus provide complementary information to that provided by eletromagnetic scatter- 
ing from the nucleus, which is sensitive only to the nuclear polar vector currents. 

Many calculational techniques have been used to determine neutrino-nuclear cross 
sections. Shell model techniques work best at lower energies where transitions to 
continuum states are not large. At intermediate energies the Continuum Random 
Phase Approximation (CRPA) is frequently used, while at still higher energies the 
Fermi gas model is thought to work well. Experimental measurements are, however, 
necessary to establish the range of validity of the different calculational techniques. 

There are more neutrino cross section measurements for 12 C than for any other 
nucleus. Three experiments, E2250 at LAMPF, the KARMEN experiment [§ at the 
ISIS facility of the Rutherford Laboratory and a liquid scintillator neutrino detec- 
tor (LSND)H ||, have measured both the exclusive reaction 12 C(is e , e~) 12 N g s and 
inclusive reaction 12 C(z/ e , e~) l2 N* to the other accessible excited states of 12 N. In 
these measurements the u e flux arises from /i + decay at rest with E u < 52.8 MeV. 
As a result of the low neutrino energy, transitions occur almost entirely to a few low 
lying states of 12 N, and over 60% of the total cross section is to the 12 N ground 
state. The cross section for producing the 12 N ground state can be predicted with 
an accuracy of ~ 5% by using model-independent form factors that can be reliably 
extracted from other measurement s[f|. All three experimental measurements of the 
12 C(z/ e , e~) 12 N g s cross section agree well with the expected value. Calculation of 
the inclusive cross section to the excited states of 12 N is model dependent and is a 
less certain procedure. The Fermi gas model (FGM) is not reliable in this instance 
because the low neutrino energy leads to momentum transfers ( Q < 100 MeV/c ), 
much smaller than the Fermi momentum ( 200 MeV/c ) in carbon. Thus extensive 
modeling of the nuclear dynamics is necessary [7|, ||, |9], [K], [11], [12[ . Final LSND 



results for this reaction^ agree with a recent shell model calculation jKJ. They are 
somewhat lower than, but consistent with, a recent CRPA calculation ||. 

Measurements also exist for two processes closely related to v e carbon scattering: 
/i~ capture on 12 C[T3 and scattering on carbon using a beam of from ir + 
decay-in-flight |14|, [15 |. For the v e carbon measurement E u pa 40 MeV, Q ~ 50 
MeV/c, and the inclusive cross section is dominated by transitions to low multipoles 
(1 + ,1~,2~). In contrast, for the carbon measurement E u pa 170 MeV, Q ~ 200 
MeV/c, and excitations occur up to 100 MeV. The fj,~ capture process, which occurs 
from the S state, is intermediate between these two processes with Q pa 90 MeV/c. 
Because these three processes occur at different energies and momentum-transfers 
they constrain different aspects of theoretical calculations. A challenging test of a 
theoretical procedure is its ability to predict all three processes. We note that the 
inclusive cross section is strongly energy and momentum transfer dependent. Thus 
the flux averaged cross section for the reaction 12 C(f At , fi~) 12 N* is approximately 200 
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Figure 1: The layout of the A6 beam stop, as it was configured for the 1993-1995 
data taking. 



times larger than the lower energy cross section for C(u e , e ) N 



The measurement fl4|, |15|| of the inclusive cross section for Civ^,^ )X several 



years ago by LSND attracted substantial interest because a CRPA calculation pre- 
dicted a cross section almost twice as large as that observed. An improved calculation 
by the same group together with an improved calculation of the neutrino energy 
spectrum and flux discussed in Section § reduced but did not eliminate the discrep- 
ancy with the measured cross section. Recent calculations using the shell model 0, [L(J 



are in better agreement with the measured cross section. In this paper we present fi- 
nal LSND results for the inclusive cross section for 12 C(v^, and for the exclusive 
reaction 12 C(V M , /i _ ) 12 A^. s .. 



2 The Neutrino Source 

The data reported here were obtained between 1994 and 1998 at the Los Alamos 
Neutron Science Center (LANSCE) primarily using neutrinos produced at the A6 
proton beam stop. As discussed below some neutrinos were also produced at upstream 
targets Al and A2. The neutrino source is described in detail elsewhere [jT~ 



In 1994 and 1995 the beam stop consisted of a 30-cm water target followed by a 
50-cm decay region, isotope production stringers and a copper beam dump as shown 
in Fig. [I]. The high-intensity 798 MeV proton beam from the linear accelerator 
generates a large pion flux from the water target. The flux of and used for 
the measurements reported here arise from the decay in flight (DIF) of n + and n~ . 
For the LANSCE proton beam and beam stop configuration n + production exceeds 
7r~ production by a factor of approximately 8 and even more for high energy pions. 
Approximately 3.4% of the n + and 5% of the 7r~ decay in flight. Upstream targets 
contributed 6% to the DIF neutrino flux. For the 1995 run, the water target was 
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removed for 32% of the 7081 C of beam. For this portion of the run the DIF z/ M flux 
was reduced approximately 50%. 

After the 1995 run the beam stop was substantially modified for accelerator pro- 
duction of tritium (APT) tests. The replacement of the water target by a close- 
packed, high-Z target resulted in reduced n + production largely due to the change 
in the neutron to proton ratio in the target. The closer packing of materials for the 
APT stop reduced the fraction of 7r =1= which decay in flight. The resulting DIF flux 
per incident proton is only one-half of that obtained with the water target. There 
were no upstream targets for almost all of the data taken with the APT target in 
1996-1998. 

The LANSCE beam dump has been used as the neutrino source for previous 
experiments |T7|, |T8| , [I9fl . These experiments primarily used neutrinos from the decay 
at rest (DAR) of stopped 7r + and A calibration experiment [ 2^] measured the 
rate of stopped /i + from a low-intensity proton beam incident on an instrumented 
beam stop. The rate of stopped fi + per incident proton was measured as a function 
of several variables and used to fine-tune a beam dump simulation program|2T|. The 
simulation program could then be used to calculate the flux for any particular beam 
dump configuration. The calibration experiment determined the DAR flux to ±7% 
for the proton energies and beam stop configurations used at LANSCE. There are 
greater uncertainties in the DIF fluxes. Uncertainties in the energy spectra of the ir^ 
which decay in flight lead to uncertainties in both the magnitudes and shapes of the 
and energy spectra. The shapes of the and energy spectra are particularly 
important for the measurement of inclusive cross sections, since these cross sections 
have a strong energy dependence. We have studied in detail the procedure of Ref. 
PTJ for calculating the DIF flux and slightly revised it to incorporate more recent 
7r + production results [E2] and to remove some slight distortions arising from use of 



finite bins in the pion production angle |23|1. The primary effect is to reduce the 
flux above 200 MeV. 

The largest uncertainties in the DIF flux arise from systematic effects in the cal- 
ibration experiment (5%), uncertainties in the ir + production cross sections used in 
the simulation (10%) and other systematic effects in the simulation (7%). For the 
1994-1995 data the upstream targets introduced a small additional uncertainty in 
the flux. For 1996-1998 data the geometry of the beam stop configuration was more 
complicated than that used in 1994-1995. The uncertainty in the DIF flux for neu- 
trinos above muon production threshold is estimated to produce an uncertainty in 
the measured cross section of 15%. This uncertainty provides the largest source of 
systematic error for the cross sections presented here. 

The LANSCE proton beam typically had a current of 800 /j,A at the beam stop 
during the 1994-1995 running period and 1000 fiA for 1996-1998. For 1994 and 
1995 the energy was approximately 770 MeV at the beam stop due to energy loss in 
upstream targets, while it was approximately 800 MeV in 1996, 1997 and 1998. 

Table [I] shows for each year the calculated flux above the threshold (123.1 MeV) 
for muon production on carbon and averaged over the LSND detector. The flux 
above threshold (113.1 MeV) for the process + p — > /i + + n is also shown. Figure 
shows the calculated and energy spectra. 
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Table 1: The and fluxes above the threshold energy for muon production and 
averaged over the LSND detector. 



Year 


flux [cm 2 


flux [cm 2 




1994 


6.04 x 10 n 


6.06 x 10 iU 


1995 


5.97 x 10 11 


6.21 x 10 10 


1996 


2.06 x 10 11 


2.91 x 10 10 


1997 


4.46 x 10 11 


5.86 x 10 10 


1998 


1.74 x 10 11 


2.32 x 10 10 


Total 


2.03 x 10 12 


2.34 x 10 11 



> 




300 



Neutrino Energy [MeV] 

Figure 2: The solid line shows the flux shape of from tt + decay- in-flight. The dashed 
line shows the flux from ix~ decay-in-flight for the same integrated proton beam. 
The muon production threshold energy for each spectrum is shown by a vertical line. 
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Figure 3: Detector enclosure and target area configuration, elevation view. 

3 The LSND Detector 

The detector is located 29.8 m downstream of the proton beam stop at an angle of 
12° to the proton beam. Figure || shows a side-view of the setup. Approximately 
2000 g/cm 2 of shielding above the detector attenuates the hadronic component of 
cosmic rays to a negligible level. The detector is also well shielded from the beam 
stop so that beam-associated neutrons are attenuated to a negligible level. Enclosing 
the detector, except on the bottom, is a highly efficient liquid scintillator veto shield 
which is essential to reduce contributions from the cosmic ray muon background to 
a low level. Reference Jl6| provides a detailed description of the detector, veto, and 
data acquisition system (DAQ) which we briefly review here. 

The detector is an approximately cylindrical tank containing 167 tons of liquid 
scintillator and viewed by 1220 uniformly spaced 8" Hamamatsu photomultiplier tubes 
(PMT) covering 25% of the surface inside the tank wall. When the deposited energy in 
the tank exceeds a threshold of approximately 4 MeV electron-equivalent energy and 
there are fewer than 4 PMT hits in the veto shield, the digitized time and pulse height 
of each of these PMTs (and of each of the 292 veto shield PMTs) are recorded. A 
veto, imposed for 15.2 /is following the firing of > 5 veto PMTs, substantially reduces 
(10 -3 ) the large number of background events arising from the decay of cosmic ray 
muons that stop in the detector. Activity in the detector or veto shield during the 
51.2 fis preceding a primary trigger is also recorded, provided there are > 17 detector 
PMT hits or > 5 veto PMT hits. This activity information is used in the analysis to 
identify events arising from muon decay. In particular, in this analysis the activity 
information is used to identify fi~ from the reaction + 12 C — > + X. For such 
events the e~ from the subsequent decay \T — > e~ + + v e provides the primary 
trigger. Activities with > 3 veto PMT hits mostly arise from cosmic ray muons and 
are rejected in the analysis. It should also be noted that the 15.2 /is veto applies only 
to the primary trigger and not to the activities preceding a valid trigger. Data after 
the primary event are recorded for 1 ms with a threshold of 21 PMTs (approximately 
0.7 MeV electron-equivalent energy). This low threshold is used for the detection 
of 2.2 MeV 7 from neutron capture on free protons. Muon events with associated 
neutrons arise from the processes u^p —>■ fi + n, v^C — > /i + nX, and v^C — > fi~nX. 

The detector operates without reference to the beam spill, but the state of the 
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beam is recorded with the event. Approximately 94% of the data is taken between 
beam spills. This allows an accurate measurement and subtraction of cosmic ray 
background surviving the event selection criteria. 

The detector scintillator consists of mineral oil (CH2) in which is dissolved a 
small concentration (0.031 g/1) of b-PBD[^4]]. This mixture allows the separation 
of Cerenkov light and scintillation light and produces about 33 photoelectrons per 
MeV of electron energy deposited in the oil. The combination of the two sources 
of light provides direction information for relativistic particles and makes particle 
identification (PID) possible. Note that the oil consists almost entirely of carbon 
and hydrogen. Isotopically the carbon is 1.1% 13 C and 98.9% 12 C. Stopping \i~ are 
captured on 12 C 8% of the time in the LSND detector. The ^ which decay are 
readily identified as muons by the presence of subsequent spatially correlated Michel 
electrons. 

The veto shield encloses the detector on all sides except the bottom. Additional 
counters were placed below the veto shield after the 1993 run to reduce cosmic ray 
background entering through the bottom support structure. More counters were 
added after the 1995 run. The main veto shield^] consists of a 15-cm layer of 
liquid scintillator in an external tank and 15 cm of lead shot in an internal tank. 
This combination of active and passive shielding tags cosmic ray muons that stop in 
the lead shot. A veto inefficiency < 10~ 5 is achieved with this detector for incident 
charged particles. 

4 Analysis Techniques 

In the analysis presented in this paper we require a fi^ followed by a delayed coin- 
cidence with a decay e . As a result of this coincidence requirement a clean beam 
excess sample of events can be obtained with relatively loose selection criteria. Fur- 
thermore, it is easy to verify that the events in this sample arise from muon decay 
since the muon lifetime and the decay electron energy spectrum are well known. 

Each event is reconstructed using the hit time and pulse height of all hit PMTs in 
the detector [TEfl. The present analysis relies on the reconstructed energy, position, and 
two PID parameters, x'tot an d « ; as described later in this section. The parameters 
Xtot an d ol are used to distinguish electron events from events arising from interactions 
of cosmic ray neutrons in the detector. Fortunately, it is possible to directly measure 
the response of the detector to electrons and neutrons in the energy range of interest 
for this analysis by using copious control data samples. We also use a Monte Carlo 
simulation, LSNDMC[^(J, to simulate events in the detector using GEANT. 

The response of the detector to electrons was determined from a large, essentially 
pure sample of electrons (and positrons) from the decay of stopped cosmic ray fj, 
in the detector. The known energy spectra for electrons from muon decay was used 
to determine the absolute energy calibration, including its small variation over the 
volume of the detector. The energy resolution was determined from the shape of the 
electron energy spectrum and was found to be 6.6% at the 52.8 MeV end-point. 

For relativistic electrons in the LSND detector approximately 65% of the photo- 
electrons arise from direct or reradiated Cerenkov light and only 35% from scintillator 
light. For muons, the threshold kinetic energy for Cerenkov radiation in the LSND 
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detector is 39 MeV. For the sample of muons analyzed in this paper only about half 
are above Cerenkov threshold and none fully relativistic. As a result, the light output 
per MeV of energy loss for the muons is significantly less than that for relativistic 
electrons. There is no calibration sample available of low-energy muons with known 
energies. Thus we rely on the Monte Carlo simulation LSNDMC for muons. We 
discuss the muon energy scale further in Sections |6| and [7] when we compare observed 
and expected energy distributions. 

There are no tracking devices in the LSND detector. Thus, event positions must 
be determined solely from the PMT information. The reconstruction process deter- 
mines an event position by minimizing a function \ r which is based on the time of 
each PMT hit corrected for the travel time of light from the assumed event position to 
the PMT[[L6|. The procedure used in several previous analyses systematically shifted 
event positions away from the center of the detector and thus effectively reduced 
the fiducial volume |3|, The reconstruction procedure has been analyzed in detail 
and an improved reconstruction procedure was developed which reduces this system- 
atic shift and provides substantially better position resolution. This procedure also 
provides results which agree well with positions obtained from the event likelihood 



procedure described in Ref ||28|| . In the analysis presented in this paper, a fiducial 
cut is imposed by requiring D > 35 cm, where D is the distance between the recon- 
structed electron position and the surface tangent to the faces of the PMTs. For the 
muon we require D > cm. 

The particle identification procedure is designed to separate particles with veloc- 
ities well above Cerenkov threshold from particles below Cerenkov threshold. The 
procedure makes use of the four parameters defined in Refs. [[I1J. Briefly, Xr and Xa 
are the quantities minimized for the determination of the event position and direction, 
Xt is the fraction of PMT hits that occur more than 12 ns after the fitted event time 
and Xtot is proportional to the product of Xr, Xa and Xt- 



Several previous LSND analyses 0, [L5|, |27j have used Xtot for particle identification. 



The distribution of Xtot for electrons, however, has a small variation with electron 
energy and with the position of the event. Therefore, in this paper, we used a modified 
variable, x'toti "with approximately a mean of zero and sigma of one, independent of 
the electron energy and positions. We also used the variable, a, which is based on 
the event likelihood procedures discussed in Ref. [2l|. It is similar to the parameter 
p discussed there, which is based on the ratio of Cerenkov to scintillator light. The a 
parameter varies from to 1 and peaks at one for electrons and at 0.3 for neutrons. 
As discussed in Refs. Q and f29j, the combination Xa = Xtot + 10(1 — at) provides 



better separation of electrons, muons, and neutrons than x'tot by itself. 

We note that a modest particle identification requirement was imposed in the 
initial data processing that created the samples analyzed here. The effect of this 
requirement is also included in the analysis. 

Cosmic ray background which remains after all selection criteria have been applied 
is well measured with the beam-off data and subtracted using the duty ratio, the ratio 
of beam-on time to beam-off time. The subtraction was performed separately for each 
year's data using the measured duty ratio for that year. The ratio averaged over the 
full data sample was 0.0632. In the present analysis the beam-off background is very 
small(< 2%) because we require a low energy muon that is well correlated in space 
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Table 2: The electron selection criteria and corresponding efficiencies. 



I Inn flflr 

quantity 


oircena 


rL/mciency 


Fiducial volume 


D > 35 cm, 


U.ooUztU.UOO 


Electron energy 


lo < hi e < ou Mev 


n oo/i -i~n m n 
U.yz4±U.UlU 


JTcLILlCie i.i-J 


Xa < 1 


n QAn-i-n m s 


In-time veto 


< 4 PMls 


0.986±0.0l0 


Future activity 


At/ > 8.8 fis 


0.99l±0.003 


Past activity 


See text 


0.858±0.0l0 


Trigger veto 


> 15.1/xs 


0.760±0.0l0 


DAQ dead time 




0.977±0.010 


Tape dead time 




0.98l±0.010 


Total 




0.467±0.032 



Table 3: The muon selection criteria and corresponding efficiencies. 



Quantity 


Criteria 


Efficiency 


Fiducial volume 


D > cm 


0.978±0.010 


Not pT Capture 




0.922±0.005 


Muon Energy (e~ equivalent) 


E < 70 MeV 


0.990±0.008 


Spatial Correlation 


Ar < 80 cm 


0.990±0.002 


H Decay Time 


0.7 < At < 9 fjs 


0.687±0.005 


Intime Veto 


< 4 PMTs 


0.988±0.010 


Total 




0.599±0.0ll 



and time with a Michel electron. 

5 Event Selection Criteria 

The analysis is designed to select the \i~ from the reaction u^+^C — > + X and the 
subsequent electron from the decay — > e~ + P e + I V In the LSND detector medium 
92% of the stopped fi~ decay and 8% are captured. The \i~ and other particles arising 
from the charge-changing neutrino interaction produce light that causes an average of 
250 PMTs to fire. The detector charge Q M , measured in photoelectrons, arises mostly 
from the \i~ but includes contributions from other particles produced in the reaction 
such as protons and 7's. 

Tables |] and ^| respectively show the selection criteria and corresponding efficien- 
cies for the electron and for the muon. For events in the decay-in-flight sample the 
event position is best determined from the reconstructed electron position rather than 
the reconstructed muon position, especially for events with low-energy muons. There- 
fore, the fiducial selection is imposed primarily on the electron. The reconstructed 
electron is required to be a distance D > 35 cm from the surface tangent to the faces 
of the PMTs. There are 3.65 x IO 30 12 C nuclei within this fiducial volume. The fidu- 
cial volume efficiency, defined to be the ratio of the number of events reconstructed 
within the fiducial volume to the actual number within this volume, was determined 
to be 0.880 ±0.055. This efficiency is less than one because there is a systematic shift 
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Figure 4: The observed energy distribution of electrons from [i~ decay for the inclusive 
sample, 12 C(z/ M , n~)X. The histogram shows the expected energy distribution of 
Michel electrons from Monte Carlo simulation. 

of reconstructed event positions away from the center of the detector, as discussed in 
Section f|. The muon is required to reconstruct only inside the region D > cm. A 
lower limit on the electron energy of 18 MeV eliminates the large background from 
12 B (3 decay created by the capture of cosmic ray on 12 C. Figure [| shows the 
observed electron energy distribution compared with the expected energy distribution 
of Michel electrons from Monte Carlo simulation. The distribution of the time, At Me , 
between the muon and electron candidates, shown in Fig. [5], agrees well with the 2.03 
fis lifetime in mineral oil. The best fit, also shown, corresponds to a lifetime of 
2.03 ± 0.05 [is. The requirement At^ e > 0.7/zs is imposed to insure that the \i and 
e are clearly separated in the trigger and in the readout of the data. The excellent 
agreement with expectations in Figs. ^ and |5| clearly shows that the events arise 
from muon decay. There is an 8% loss of events due to \i~ capture in the detector 
medium. Figure |] shows the spatial separation Ar between the reconstructed muon 
and electron positions. A loose requirement, Ar < 0.8 m, is imposed to minimize the 
background from accidental fi, e correlations while retaining high acceptance. 

Many of the selection criteria are designed to reduce the cosmic ray background, 
especially that arising from the decay of cosmic ray muons which stop in the detector. 
Both the muon and the electron candidates are required to have fewer than 4 PMT 
hits in the veto and no bottom counter coincidence during the 500 ns event window. 
The detector PMT faces are 25 cm inside the tank and thus stopping cosmic ray muons 
must traverse at least 60 cm of oil to reach the fiducial volume. As a result, these 
muons typically produce a large detector signal. The requirement E <70 MeV, where 
E is the electron equivalent energy of the muon, eliminates most such background 
events with almost no loss of acceptance for muons arising from neutrino interactions. 
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Figure 5: The distribution of the time difference, At Me , between the \x~ and the decay 
e~ in the inclusive sample, 12 C(z/ M , /j,~)X . The best fit (solid line) curve corresponds 
to a lifetime of 2.03 ± 0.05 /is. 




Figure 6: The distribution of the distance between the reconstructed positions of the 
/i~ and e~ in the beam-excess inclusive sample, 12 C(z/ M , fJ,~)X. The histogram is the 
prediction from the Monte Carlo simulation, normalized to the data. 
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Figure 7: The distribution of the particle identification parameter Xa of electrons 
from fi~ decay for the inclusive sample 12 C(z/ A1 , /i _ )X. The histogram shows the Xa 
distribution of Michel electrons. 



Muons which are misidentified as electrons are removed by requiring that there 
be no future activity consistent with a Michel electron. Any electron candidate with 
future activity with fewer than 4 veto PMT hits and more than 50 detector PMT hits 
within 8.8 /is is rejected. 

Frequently, in addition to the candidate muon which satisfies the criteria in Table 
there are one or more other activities prior to the electron. If an activity is due to 
a stopping muon, that muon could be the parent of the observed electron. Therefore 
an event is rejected if, in the 35 /j,s interval prior to the electron, there is an activity 
(other than the muon) with Q > 3000 pe or an activity with > 4 PMT hits in the 
veto and > 100 PMT hits in the detector. We also reject any event with a past 
activity within 51 fis with > 5 veto PMT hits and > 500 detector PMT hits. A 
cut is performed during initial data processing on past activities that are spatially 
correlated with the primary event, within 30 fxs of the primary event, and have > 4 
veto PMT hits. 

The acceptances for the past activity and in-time veto cuts are obtained by ap- 
plying these cuts to a large sample of events triggered with the laser used for detector 
calibration. These laser events are spread uniformly through the run and thus average 
over the small variation in run conditions. 

Only a loose particle ID requirement, Xa < 4.0, was imposed on the electron 
and none on the muon. A sample of Michel electrons (electrons from the decay of 
stopped /r* 1 ) was analyzed to obtain the acceptance of electrons for the Xa particle 
identification cut, as shown in Fig. [7|. 

The Monte Carlo simulation LSNDMC was used to obtain the PMT hit distribu- 
tions expected from the various processes that contribute to the inclusive sample and 
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PMT Hits 



Figure 8: The observed PMT hit distribution for the decay-in-flight sample (including 
Vfj, + C — > fi~ + X, u^C — > fi + X and u^p — > /i + n). The histogram is the prediction 
from the Monte Carlo simulation, normalized to the data. 



to the exclusive sample with an identified (3 decay. Figure [8] compares the observed 
and expected distribution of PMT hits for muons from inclusive events. There is ex- 
cellent agreement, and thus we expect that the simulation provides a reliable estimate 
of the fraction of events lost because they are below the activity threshold of 18 PMT 
hits (roughly 4 MeV). For the inclusive sample (exclusive sample) we find that only 
1.2%(2.4%) of the events have fewer than 18 PMT hits. 



6 The Transition to the 12 N Ground State 

For analysis of the exclusive process 12 C(z/ /U , fi~) 12 N g , s , we also require detection of 
the e + from the (3 decay of 12 N g s . Therefore, for these events three particles are 
detected: the muon, the decay electron and the positron from the (3 decay of 12 N, 



g.s. 



Table [| gives the selection criteria and efficiencies for the 12 N (3 decay positron. These 
are the same criteria used previously in an analysis of a much larger sample from the 
analogous process 12 C(z/ e , e~) 12 N gs || |J. The f3 decay has a mean lifetime of 15.9 ms 
and maximum positron kinetic energy of 16.3 MeV|3(|. Figure ^ shows the observed (3 



decay time distribution compared with the expected 15.9 ms lifetime. Figure |10| shows 
the distance between the reconstructed electron and positron positions for the beam- 
excess sample. A cut was applied at 70 cm, resulting in an acceptance of (99.2±0.8)%. 
The positron is required to be spatially correlated with the electron rather than the 
muon because the position of the electron in general is better determined. Following 
a muon produced by a neutrino interaction, an uncorrelated particle, such as the 
positron from 12 B [3 decay, will occasionally satisfy all the positron criteria including 
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Table 4: Criteria to select e + from N g . s . beta decay and corresponding efficiencies for 
the reaction 12 C(^, /i~) 12 iV 9 . s .. 



Quantity 



Criteria 



Efficiency 



(3 decay time 
Spatial correlation 
PMT threshold 

Positron energy 
Fiducial volume 
Trigger veto 
In-time veto 
DAQ dead time 



52 /is< t < 60 ms 
Ar < 0.7 m 

> 100 for 1994, 

> 75 after 1994 
E/3 < 18 MeV 

D>0cm 
> 15.1 /js 
< 4 PMTs 



0.974±0.002 
0.992±0.008 
0.856±0.011 

0.999±0.001 
0.986±0.010 
0.760±0.010 
0.988±0.010 
0.977±0.010 



Total 



0.598±0.016 




Figure 9: The distribution of time differences between the electrons and (3 in the 
exclusive sample of 12 C(z/ At , /j,~) 12 N g . s . is compared with the expected (3 lifetime. The 
dotted line shows the calculated accidental contribution. 
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Figure 10: The distribution of the distances between the electrons and (3 for beam- 
excess events in the exclusive sample of 12 C(^ M , n~) 12 N g , s ,. The dashed line shows the 
estimated accidental contribution. The solid line shows the expected shape (including 
the accidental contribution) from the Monte Carlo simulation, normalized to data. 

the requirements of time (60 ms) and spatial (70 cm) correlation with the electron. 
The probability of such an accidental coincidence can be precisely measured from the 
Michel electron sample. The background from this source is also shown in Figs. [| 



and [10]. The efficiency of 76.0% caused by the 15.2 fxs veto and the trigger dead time 
of 2.3% are the same as for the electron. Positrons with four or more in-time veto hits 
or any bottom veto coincidence are rejected. The Monte Carlo simulation was used 
to generate the expected distribution for the positron energy. There was a trigger 
requirement of 100 PMT hits for 1994 and 75 PMT hits for 1995-1998. The positron 



was required to have an energy less than 18 MeV. Figure [LI] compares the observed 
and expected positron energy distributions. Figure |T2] compares the observed and 
expected distributions of muon decay time. 

Excited states of 12 N decay by prompt proton emission and thus do not feed 
down to the 12 N ground state or contribute to the delayed coincidence rate. The 
form factors required to calculate the cross section are well known from a variety 
of previous measurements ||. This cross section and the known flux are used to 
obtain the expected muon kinetic energy spectrum which is compared with the data 
in Fig. [TS|. 

As stated in Section |4] the energy calibration for muons (the conversion from 
photoelectrons to MeV) is obtained from the Monte Carlo simulation LSNDMC. For 
this ground state reaction, the expected muon energy distribution should be very 
reliable. Thus the agreement seen in Fig. [13] provides confirmation for the muon 
energy calibration within the limited statistics. 
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Ep [MeV] 

Figure 11: The distribution of (3 energy from the exclusive sample of 
12 C(z/ At , fi~) 12 N g , s ,. The dashed line shows the estimated accidental contribution. The 
solid line shows the expected shape (including the accidental contribution) from the 
Monte Carlo simulation, normalized to data. 




Figure 12: The distribution of muon decay time from the exclusive sample of 
12 C(z/ M , /i~) 12 N g _ s _. The curve shows the distribution from the expected muon life- 
time in the oil. 
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[MeV] 



Figure 13: The observed and expected (solid line) \x kinetic energy distribution for 
beam excess-events in the 12 C(f M , n~) l2 N g . s . sample. 



There are two sources of background. The largest arises from the accidental coin- 
cidence of a positron candidate with an event from the inclusive sample of neutrino- 
induced muons. The probability of an uncorrelated particle satisfying all the positron 
criteria, including the requirements of time (60 ms) and spatial correlation (70 cm) 
with the electron, can be precisely measured from a large Michel electron sample. 
The second background arises from the process 12 C(v^, /i + ) 12 -B g . s ., where we detect 
the e~ from the f3 decay of the 12 B ground state[|JT|. This background is small pri- 
marily because the flux of high-energy is approximately a factor of 10 lower than 
the corresponding flux and because the 12 B g lifetime is longer than the l2 N g . s , 
lifetime. 

Table |5| shows the number of beam excess events, the number of background 
events, the efficiency, the neutrino flux for E v > 123.1 MeV, and the cross section 
averaged over the flux. The efficiency shown includes the efficiency for muons to 
have more than 17 PMT hits as well as the electron, muon and beta efficiencies given 
in Table Table |3| and Table f|, respectively. The flux-averaged cross section is 
(a) = (5.6 ± 0.8 ±1.0) x 10~ 41 cm 2 , where the first error is statistical and the second 
systematic. The two dominant sources of systematic error are the neutrino flux (15%) 
discussed in Section |2] and the effective fiducial volume (6%) discussed in Section |. 

For this reaction to the l2 N ground state it is also straightforward to measure the 
energy dependence of the cross section. The recoil energy of the nucleus is small and 
to a very good approximation, E u = m^c 2 + T M + 17.7 MeV, where m M is the muon 
mass, T M is the kinetic energy of the muon, and 17.7 MeV arises from the Q value of 
the reaction and the nuclear recoil. Figure 111 compares the measured cross section 



as a function of E v with three theoretical calculations obtained from Ref. 31 . The 
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Table 5: Beam-excess events, background, efficiency, neutrino flux and flux-averaged 
cross section for the exclusive reaction 12 C(u^, fi~) 12 N, 



g.s. 



Corrected beam excess events 
v, ± 12 C - fi + + 12 B g . s . 
accidental e + background 



77.8±8.9 
2.7±0.5 
8.2±0.8 



+ 12 N 



12 c 

Efficiency 

z/ M flux (E u > 123.1 MeV) 
(a) measured 



66.9±9.0 
16.3 ± 1.2% 
2.03xl0 12 cm" 2 
(5.6 ±0.8 ±1.0) x 10" 41 



cm 



(a) theory 



P3 



Engel et al. 
Kolbe et al. § 
Volpe et al. @ 
Hayes and Towner |T0 



6.4 x 10" 41 cm 2 



41 cm 2 



7.0 x 10 

6.5 x 10 

5.6 x 10" 41 cm 2 



41 cm 2 




150 200 250 

Neutrino Energy(MeV) 



Figure 14: The measured cross section for the process 12 C(z/ A1 ,/i" 
with three theoretical calculations obtained from Ref. [0111. 



\12 



N gs compared 
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Table 6: Beam-excess events, background, efficiency, neutrino flux and flux-averaged 
cross section for the inclusive reaction 12 C(u fl , ^) 12 X. 



\,c\WP>oi'P>f\ hpain pvppcici pa^pti^c; 


2464+5D 

l—t 1 rill 


v^+p ^ [i + + n 


217±35 


P„ + 12 C -> /i + + X 


71±35 


^ + 13 C -> M~ + X 


24±12 


^ + 12 c -> /x- + X 


2152±56 


Efficiency 

^ flux (£„ > 123.1 MeV) 


(27.7 ± 1.9)% 
2.03xl0 12 cm" 2 


(a) measured 


(10.6 ± 0.3 ± 1.8) x 10" 40 cm 2 


Theory 


Kolbe et al. J§ 


17.5 x 10~ 40 cm 2 


Volpe et al. § 


15.2 x 10~ 40 cm 2 


Hayes and Towner [10 


13.8 x 10~ 40 cm 2 



agreement is reasonable within the limited statistics. 

There is little disagreement ||31|| in the predicted cross section for this exclusive 
process, as it is fixed by measured values of closely related electroweak transition 
probabilities. The differences that exist among various calculations result from differ- 
ent models for the dependence of various elements of the transition probability on the 
momentum transfer. Thus, as shown in Fig. [14], the differences among shell model 
approaches, an RPA calculation, and an "elementary particle" model all agree for E v 
up to 160 MeV and differ only by about 10% at 250 MeV. 

Table || also shows three more recent calculations || [9|, [10J of the cross section 
for this exclusive process. The focus of these papers, however, was on obtaining a 
satisfactory overall description of neutrino reactions on carbon, including the inclusive 
cross section and of pr capture. 



7 The Inclusive Reaction 

Most of the inclusive beam-excess events arise from the reaction 12 C(v^, p~)X , but 
approximately 10% are due to other sources. Table [] shows the number of beam- 
excess events, the calculated backgrounds, the efficiency, flux, and the flux-averaged 
cross section for this process. The efficiency shown includes the efficiency for muons 
to have more than 17 PMT hits as well as the electron and muon efficiencies given in 
Tables ^ and |||. The backgrounds arising from the component of the decay- in-flight 
beam are small, primarily because the high-energy flux is approximately a factor 
of 10 lower than the corresponding flux. The largest background arises from the 
process + p — > /i + + n. The cross section is well known and the uncertainty in 
this process is mainly due to the 15% uncertainty in the flux. A much smaller but 
less well understood background arises from the process 12 C(P M , p, + )X. Plausibly, 
as observed for the process 12 C(z/ A1 , p,~)X, the cross section might be expected to 
be approximately 60% of that given by a recent CRPA calculation^!. We use this 
reduced cross section in calculating this background but assign a large error to reflect 
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the uncertainty in the cross section. An even smaller background arises from the 1.1% 
13 C component of the scintillator. For the process 13 C(z/ M , fi~)X we use a Fermi gas 
model calculation and assign a 50% uncertainty. 

The measured flux-averaged cross section for the inclusive reaction 12 C(z/ M , n~)X 
is (a) = (10.6 ± 0.3 ± 1.8) x 10~ 40 cm 2 , where the first error is statistical and the 
second systematic. The mean energy of the neutrino flux above threshold is 156 MeV. 

Table |] shows for comparison results of three recent theoretical calculations of 



the cross section R || 10]. The shell model calculation of Hayes and Towner [1C] 



provides the best agreement with our measurement but all three predict a higher cross 
section than observed. These calculations used the energy distribution calculated 
by LSND for the 1994 beam stop. The shape of the calculated energy distribution, 
however, shows only small yearly variations that would produce shifts of a few percent 
in the expected cross section for the different years. Averaged over all the data the 
calculated shape is the same as that calculated for 1994. 

The systematic error is due almost entirely to the uncertainty in the flux. 
To determine this value, the inputs to the neutrino beam Monte Carlo program 
were varied within their estimated uncertainties. The resulting variation in both the 
magnitude and the shape of the flux above muon production threshold results in 
a 15% uncertainty in the inclusive cross section. 

The beam dump configuration was substantially modified after the 1995 run as 
discussed in Section |2|. The two beam dump geometries are very different. Also the 
flux for 1994-1995 depends primarily on 7r + production from water while the flux for 
1996-1998 arises primarily from 7r + production on high-Z materials. The resulting 

DIF flux per incident proton with the modified configuration was only one-half of 
that obtained with the water target used in 1994-1995. Thus a comparison of the 
cross sections measured with the two beam dump configurations provides a check 
on the systematics of the flux simulation. The cross section measured with the 
APT beam dump in 1996-1998, (9.5 ± 0A(stat.)) x 10" 40 cm 2 , was lower than that 
measured with the water target setup in 1994-1995, (11.4 ± 0A(stat.)) x 10 -40 cm 2 , 
by a little more than one standard deviation. The statistical uncertainties of these 
measurements are small compared to the systematic uncertainties. The systematic 
errors for the two beam dump configurations are partially correlated but each has 
an uncorrelated systematic error of about 10%. The cross sections for the two beam 
dump configurations differ by 1.9 x 10 -40 cm 2 compared to an expected uncertainty 
of 1.5 x 10 -40 cm 2 . Thus results for the two beam dump configurations are consistent 
within errors. 



The spatial distribution of the beam-excess electrons is shown in Fig. [L5[ There 
is a clear enhancement of events at high x and high y due to the variation of the 
flux over the detector. The good agreement with expectation shows that this spatial 
distribution is well modeled by the beam simulation program. 

For the reaction 12 C(v^, n~)X , the detector charge Q^, measured in photoelec- 
trons, arises mostly from the pT but includes contributions from other particles in 
the reaction such as protons and 7's. The muon kinetic energy distribution obtained 
from the distribution is shown in Fig. |16[ According to the CRPA calculation 
32] , 84.5% (20.4%) of the inclusive 12 C(z/ M , p~)X events have proton (7) emission 



from the decay of the N excited states. The average energies of these protons and 7s 
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Figure 15: The spatial distribution of the electron for beam-excess events compared 
with expectation (solid line) from the inclusive reaction 12 C(z/ M , fJ,~)X. 




Figure 16: The observed and expected (histogram) distributions of the muon kinetic 
energy, E^, for the inclusive decay-in-flight sample. The expected distribution has 
been normalized to the data. 
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Table 7: The expected and observed numbers of events with associated neutrons and 
the calculated background from reactions. 





Events from 


Fraction 


Events with 


Source 


Table |6| 


with neutron 


neutron 


Beam Excess 


2464 


(11.6 ± 1.1)% 


286±27 




214 


100% 


214±35 


v^C — > n + nX 


64 


79% 


51±26 


v^C -> \i~X 


2181 


1.0 ±2.4% 


21±52 



are estimated to be 9 and 4 MeV. We used the calculation of Ref. |3j to determine 
proton and 7 energy distributions and LSNDMC to determine the number of pho- 
toelectrons produced. The averaged contribution of these particles to the measured 
electron equivalent energy of the muons is estimated to be 2.2 MeV. Protons produce 
less scintillation light than electrons due to saturation effects. The uncertainty in the 
saturation effect is the primary source of uncertainty in the muon and proton energy 
determination. The average contribution to from particles other than the muon is 
approximately 12%. The expected muon energy spectrum in Fig. [16] is obtained from 
a recent CRPA calculation and includes the contribution from particles other than 
muon. There is reasonable agreement. However, given the uncertainties in the shape 
of the energy spectrum, in the modeling of the energy from nuclear breakup and 
in the muon and proton energy calibration, we do not try to extract any information 
on the energy dependence of the cross section for the reaction 12 C(z/^, fJ>~)X. 

Further information on the inclusive sample can be obtained by measuring the 
fraction of the events with an associated neutron. The presence of a neutron is 
established by detection of the 7 ray from the neutron's capture on a proton in the 
detector via the reaction n + p — > d + 7. A detailed discussion on the procedure used 
to detect these 7's can be found in Ref. ||33| . The distribution of the likelihood ratio 
R for correlated 7's from neutron capture is very different from that for uncorrelated 
(accidental) 7's. The measured R distribution for the inclusive sample, shown in Fig. 
|i~7| , was fit to a mixture of these two possible 7 sources to determine the fraction of 
events with associated neutrons. The best fit, also shown in the figure, corresponds 
to a fraction of events with a neutron of (11.6 ± 1.1)%, where the error includes 
systematic uncertainties. 

The two largest backgrounds shown in Table |5] arise from the component of the 
beam, and almost all of these events should have an associated neutron. In contrast, 
most of the events arising from interactions will not have an associated neutron. A 
CRPA calculation predicts that 79% of the events from the reaction 12 C(z^, n + )X will 
have an associated neutron compared to only 5.9% for the reaction 12 C(u^, 

Table shows the measured component with an associated neutron for the beam 
excess sample and the calculated backgrounds from interactions. The resulting 
number for the carbon sample (1.0 ± 2.4%) is lower than the CRPA prediction of 
5.9%. The CRPA prediction is for the spectrum used in our earlier analysis fl5|l , 
which is slightly harder than that used in the present analysis as discussed in Section 
0. The observed number of events with neutrons also rules out a flux much bigger 
than that calculated by the beam Monte Carlo simulation. 
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Figure 17: The observed distribution of the 7 likelihood ratio R for the inclusive 
decay-in-flight sample. Shown for comparison are the best fit (solid line) combination 
of the correlated distribution and uncorrelated distribution to the data. The best fit 
has a (11.6±1.1)% correlated component. The shaded region shows the uncorrelated 
contribution. 



8 Conclusions 

The exclusive process 12 C(z^, yU~) 12 iV 5 . s . has been measured with a clean sample of 
66.9 ± 9.1 events for which the fi~ , the decay e~, and the e + from the (3 decay of the 
12 N g Si are detected. For this process the theoretical cross section calculations are very 
reliable. The flux averaged cross section is measured to be (5.6 ± 0.8 ±1.0) x 10 -41 
cm 2 in good agreement with theoretical expectations. From comparison of this cross 
section with the cross section for the inclusive process 12 C(^ At , n )X we obtain a 
flux-averaged branching ratio of (5.3 ± 0.8 ± 0.5)%. 

The inclusive process 12 C{p^^~)X has also been measured. There are model- 
dependent uncertainties in the theoretical cross section calculation for this process 
that are not present for the 12 N g , s , cross section. The measured flux-averaged cross 
section, (10.6±0.3±1.8) x 10~ 40 cm 2 , is smaller than theoretical expectations [3, |9], |T0]j . 
It is in better agreement with the shell model calculation of Hayes and Towner [|l(J 
than with the CRPA calculation of Kolbe et al. M. The mean energy of the neutrino 
flux above threshold is 156 MeV. The measured distribution of the muon energy 
(including contributions from other particles such as protons and 7's) agrees within 
errors with the CRPA calculation. The fraction of events with associated neutrons 
was measured to be (1.0 ± 2.4)%. This is lower than the CRPA calculation of 5.9%. 
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